Skyrmion-like excitations in dynamical lattices 
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We construct discrete analogs of Skyrmions in nonlinear dynamical lattices. The Skyrmion is 
built as a vortex soliton of a complex field, coupled to a dark radial soliton of a real field. Adjusting 
the Skyrmion ansatz to the lattice setting allows us to construct a baby-Skyrmion in two dimensions 
(2D) and extend it into the 3D case (ID counterparts of the Skyrmions are also found). Stability 
limits for these patterns are obtained analytically and verified numerically. The dynamics of unstable 
discrete Skyrmions is explored, and their stabilization by external potentials is discussed. 



Introduction. Recently, studies of intrinsic localized 
modes (ILMs) in nonlinear lattice systems have drawn 
much attention JJ due to their relevance to various phys- 
ical problems including, inter alia, optical waveguide ar- 
rays 3 , photonic crystals p] , Bose-Einstein condensates 
(BECs) trapped in deep optical lattices (OLs) Q, and 
Josephson-junction ladders [3- 

A wide variety of self-supporting ILM type of exci- 
tations have been predicted theoretically and observed 
experimentally. Prominent examples are bright and dark 
optical discrete solitons 0, Q in AlGaAs waveguide ar- 
rays H, multi-dimensional solitons in photonic lattices 
H, Igl lldl. discrete vortex solitons E3j IT^- l attice 
dipole solitons [Ij|, multi-component solitons [l^, soli- 
ton trains [la|, necklace solitons [l3|, and so on. Paral- 
lel to these achievements in optics, a remarkable recent 
result was the creation of gap solitons in BECs loaded 
in OLs [13 ■ In isict, BECs offer an interesting imple- 
mentation of nonlinear dynamics in confined lattices, due 
to the fact that experiments are always run in an ex- 
ternal trap [3|- These developments raise the question 
whether counterparts of more complex structures, that 
were originally proposed in continuum media within field- 
theoretical contexts, can be predicted and observed in 
dynamical lattices. 

Among the most fundamental objects of this type are 
the three-dimensional (3D) Skyrmions, initially proposed 
to explain topologically the origin of the baryon quantum 
number [l9j. On the other hand, their two-dimensional 
(2D) version, so-called "baby-Skyrmions" , were used to 
model bubble generation in condensed-matter systems 
in the presence of an external magnetic field. In par- 
ticular, baby-Skyrmions seem to have a central role in 
the dissapearance of antiferromagnetism and the onset of 
high-Tc superconductivity ,2Q|, as well as in the ground 
state properties of the quantum Hall ferromagnets [21| . 
Recently, stable Skyrmions have also been predicted in 
BECs ^j. On the other hand, as concerns discrete sys- 
tems, stable discrete Skyrmions on lattices were consid- 
ered in the framework of the Heisenberg model for mag- 
netism [23 , as well as in the 2D Hubbard model |2J] . Ad- 
ditionally, Skyrme lattices are also very interesting due 
to their ability to describe electron spin textures in quan- 



tum Hall systems [23 ■ Importantly, discrete Skyrmions 
play also an essential role for the quantization of the orig- 
inal Skyrme model |26l | , the continuum version of which 
is non-renormalizable. 

In the present work, we discuss the discretization of 
Skyrmions from a different point of view. Putting the 
continuity aspect on second priority, we develop discrete 
solitons with topological properties equivalent to those 
of Skyrmions in the continuum field theory. We show 
that the proposed Skyrmion-like structures (for simplic- 
ity called "Skyrmions" hereafter) can be stabilized on the 
dynamical lattice under consideration, namely a nonlin- 
ear Schrodinger lattice. 

We aim to construct discrete Skyrmion states in a 
paradigmatic nonlinear-lattice model, viz., the discrete 
nonlinear Schrodinger (DNLS) equation. DNLS is a uni- 
versal envelope wave equation for a variety of Hamilto- 
nian systems (such as, e.g., nonlinear Klein-Gordon lat- 
tice models), and, moreover, is a direct model for BECs 
trapped in strong OLs jj| and crystals built of microres- 
onators [23 ; in addition, its 2D version models optical 
waveguide arrays |^ |y]. As Skyrmions necessarily in- 
volve (through the so-called hedgehog ansatz 28, 29]) 
three scalar fields, and are characterized by two indepen- 
dent topological charges, their description requires two 
complex field variables. Thus, in this work we study a 
two-component DNLS equation. This model is directly 
relevant to waveguide arrays jlSi] . when the light has dif- 
ferent polarizations or frequencies, and to binary BECs 
composed as spin state mixtures of the same isotope 30] . 

The paper is structured as follows: In the next section, 
we present the model and some analytical considerations. 
Then, we report numerical results for ID, 2D, and 3D 
lattices, and, finally, we summarize our findings. 

Model and Analytics. We consider a vector DNLS 
equation for <j) = (0i, 4>2)'^ on the cubic/square lattice. 
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where overdot stands for derivative in time (or, in the 
optical- waveguide model, in the propagation distance), C 
is a coupling constant, the parameter f3 characterizes the 
intra-species interaction, and the £)-dimensional discrete 



Laplacianis A^^^cj)^ = X]m=<n> 0m — 2I?0n, with < n> 
standing for the nearest-neighbor shell of n [the latter is 
a vector index in the D-dimensional case — e.g., n = 
(rii,n2,n3) in 3D]. Stationary solutions are looked for as 
(/)n = exp(— jAi)(un, Wn)"^, whcrc A is the frequency (or 
chemical potential in the context of BECs), and Un,Vn 
obey the following equations 

-CA7/„ + (|m„P+/J|^„|2)„„, (2) 

-CA«„ + (|t;„|2+/3|u„ni;„. (3) 
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Note that the nonlinear-interaction matrix G in Eq. 
(Q implies self-repulsion of the lattice field, while the 
model with self-attraction, i.e., G replaced by — G, can 
be transformed into the present form by the known 
staggering transformation; for instance, («„, ^n) = 
(-l)"i+"2+"3 (u„,w„) in the 3D case. 

We first construct solutions with the desired structures 
in the so-called anti-continuum (AC) limit of C ~ 0, 
and then continue solution branches to C > 0, by 
means of a fixed-point iteration converging to the rele- 
vant solutions. Once these have been obtained, linear- 
stability analysis is performed for a perturbed solution, 
0(P"*) = (0„ + a„e^* + b„e^**) e-'^\ where an and b„ 
constitute an eigenmode of infinitesimal perturbations, 
and A is the corresponding eigenvalue. The stationary 
solution is unstable if at least one pair of A has nonzero 
real part. 

We start with the 2D case, where the baby-Skyrmion 
can be found in the continuum model, using a "hedge- 
hog" ansatz (in the polar coordinates) [23, 123 , namely 

* = [sin(/(r))cos(A:0), sin(/(r)) sin(A;0), cos(/(r))]^ , 

(4) 
subject to the boundary conditions limr^o/(^) = Mn, 
and lim^^oo f{r) — Ntt, with integer k, M and N. The 
winding number of the corresponding state is [(—1)^ — 
(— l)*^]fc/2. We will try to construct a counterpart of 
ansatz Q) on the 2D lattice. In particular, combin- 
ing the first two components of ansatz Q as ^i -I- i^2, 
we obtain a complex field with vorticity k. To consider 
the fundamental Skyrmions, we set fc = 1 and M = 0, 
A^ = 1. Accordingly, the complex field is a localized 
discrete vortex, while the remaining field may be real, 
taking the form of a dark soliton in the (quasi-)radial di- 
rection on the lattice. For C = 0, one can construct a 
discrete vortex (alias vortex cross) by assigning the com- 
plex fields phases 0, 7r/2, tt, 37r/2 and modulus \/A at four 
sites, (TO,n) = (1,0), (0, l)j_(-l, 0) and (0,-1), which 
surround the central point U^. The radial dark soliton 
in the real field is achieved, in the same AC limit, by set- 
ting uo,o = 1 Bind Vm.n = — 1 at all other sites, with the 
exception of the above-mentioned four sites surrounding 
the origin, where Vm,n = 0. Notice that this corresponds 
to a hedgehog ansatz with /(O) = 0, /(I) = tt/2 and 
/(r > 1) = TT. Of course, the size of the vortex structure 
may be made larger in the AC limit, but this adversely 
affects the stability of the vortex [31| , therefore we do not 
examine such cases here. 



We note that the above 2D structure suggests a ID 
analog of the Skyrmion, essentially as a cross-section of 
the 2D profile. In particular, in the AC limit, the ID 
configuration becomes u„ = vA(<5ti,i — ^„,_i), vq — vA, 
vi = V-i = and Vn = — VA for \n\ > 1. This ID struc- 
ture essentially consists of a twisted localized mode |23| 
in the one field, coupled to a pair of discrete dark solitons 
in the other. A variety of generalizations to the 3D case, 
resulting from Skyrmion's rotation, are possible, but we 
limit our considerations to the simplest setting where Un 
is a planar excitation (as above) with the well-defined 
vorticity around the third dimension, while wo,o,o ~ v^, 
the field vanishes at the immediate neighbors of (0, 0, 0) 
and has a value of — vA elsewhere. 

To examine instabilities relevant to the discrete 
Skyrmion, we first examine the dispersion relation for 
the equations linearized around such a solution (i.e., the 
continuous spectrum). While the full spectrum also con- 
tains isolated eigenvalues, it is its continuous component 
which is primarily responsible for instabilities (see be- 
low). Thus, we perturb the asymptotic lattice field, far 
from the Skyrmion's center, as follows: 

02 = -l + eC'e'("*+'=")-HeDe-*('^**+'="l (6) 

Using Eqs. (0)-® in Eq. ^ yields two excitation 
branches, associated, respectively, with the zero and non- 
vanishing backgrounds of the first and second fields. 



= ± [A - /3 - ACD sin2(fc/2)] , 



uj^±JA + 4CD sin^(fc/2) - A2, 



(7) 
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which restricts the relevant spectral 

ban ds to ±[A - /3 - iGD, A - /3] and 
[- v/(A-h4Ci:»)2-A2, y/{AT4CDy^n^] . These 

bands have opposite Krein signatures J33, hence their 
collision, occuring with the increase of C, at 



C(f) = (/3-A)V[8i?(2A-/3)], 



(9) 



will generate complex eigenvalues, i.e., instability. Thus, 
the discrete Skyrmions may only be stable in the interval 
of < C < Ccr ; in particular, for Ccr > one requires 
2A > /3. It is interesting to mention that Ccr -^ oo (the 
continuum limit) as A — P/2 — > -1-0; however, analysis of 
this special case is beyond the scope of this paper. 

Numerical Results. We first display numerical find- 
ings for the ID case in Fig. ^ Shown are norms of the 



solution, A^i = X;„ |wnP and A^2 = E„(A 



), and 



most unstable eigenmodes and eigenvalues (computed on 
a lattice with 400 sites). In this case, Eq. © predicts 
Ccr ~ 0.102083, which coincides with the numerical 
finding (0.102 ± 0.0005). Examples of stable and un- 
stable ID Skyrmions are included too, for C = 0.05 and 
G = 0.15. As expected (see above), the destabilization 
occurs, indeed, via the collision of two continuous bands 
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FIG. 1: (Color online) ID discrete analog of the Skyrmion. 
Top and bottom panels in the left column show, respectively, 
norms A'^i , A'2 (defined in the text) and two most unstable 
eigenmodes vs. coupling C [the vertical dashed line is the in- 
stability onset as predicted by Eq. Q]- Upper and bottom 
panels in the middle and right columns display, respectively, 
the solution for C — 0.05 and C = 0.15, and the correspond- 
ing spectral planes (Ar, Xi) of the eigenvalue A = Ar + iXi- 
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FIG. 2: (Color online) Same diagnostics as in Fig. 0for dis- 
crete 2D (baby-) Skyrmions. Left and right paired panels 
pertain, respectively, to C = 0.025 and C — 0.075. Solution 
profiles are shown on the right, with the top, middle, and bot- 
tom rows displaying countours of real and imaginary parts of 
the first (complex) field, and the second (real) field. 



of eigenvalues at the critical point. It has been checked 
that the numerical results, displayed here ioi A — 2 — 8/3, 
adequately represent a large area in the parameter space. 

Results for the 2D case are presented in Fig. |2] We 
again observe excellent agreement of the theoretical pre- 

diction for the instability onset, at C = CcV — 
0.051042. Typical solution examples are shown for C = 
0.025 (stable) and C = 0.075 (unstable). 

An example of one of the possible (as mentioned above) 
3D generalizations of the 2D discrete baby-Skyrmion is 
shown in Fig. O In this case, the complex field is ar- 
ranged as a 3D soliton carrying a vortex in the horizontal 
plane (cf. Ref. i34\), while the real component features 
a 3D radial dark soliton. A stable 3D discrete Skyrmion 

is shown for C = 0.01 (in this case, C^? = 0.034028). 

The next step is to simulate the evolution of unsta- 
ble solutions. As, beyond the primary instability thresh- 
old @, a cascade of secondary instabilities is produced 
by additional collisions between the continuous spectral 
bands, one may expect that the corresponding multi- 
tude of unstable eigenmodes leads to a kind of lattice 
turbulence, especially because the unstable eigenmodes 





FIG. 3: (Color online) A 3D discrete Skyrmion for C = 0.05. 
The left and right panels show, respectively, contours of the 
first (complex) field, at Ke{(j)i) — ±1 (blue/red) and Im((;6i) = 
±1 (green/yellow), and of the second (real) field, at (j)2 = 
(—1,0, 1) (red, green and blue, respectively). 
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FIG. 4: (Color online) Instability development in one- and 
two-dimensional discrete Skyrmions. The top, middle, and 
bottom left panels show, respectively, space-time contours of 
the absolute value of the first (complex) field and of the second 
(real) one, and the spatial distribution of the fields at i = 1100 
in the ID case. The top right panels display the contours of 
the square modulus of the two fields for the 2D case, at the 
instability onset, t — 220, while the panels beneath them 
show the evolution of the fields at the central site and its four 
neighbors, for the same case. 



are delocalized, as they belong to the continuous spec- 
trum. Indeed, this is what we observe in direct sim- 
ulations, as shown in Fig. 0] for C = 0.149 > Ccr and 
C = 0.099 > Ci? in the ID and 2D cases. In the ID case 
(the left part of the figure), the weakly unstable config- 
uration remains undisturbed for a long time, but even- 
tually, around t w 650, the instability generates spatial 
chaos in the real (second) field component, and breath- 
ing in its complex counterpart. This dynamics persists 
for long evolution times. Transition to chaotic behavior 
is also observed (in the right part of Fig. 0J in the 2D 
case, after the onset of the instability around t = 220. 

Finally, we also considered the influence of external po- 
tentials, which is necessary in the application to BECs, 
that are typically confined by a parabolic trap. In the 
ID case, the latter amounts to adding a term {fl^ / 2)n'^ (j) 
to Eq. Q (in the 2D and 3D cases, the trap produces 
qualitatively similar effects). As seen in Fig.jSJ the profile 
of the V field now has a finite size, as per the Thomas- 
Fermi approximation 4J . The main novel feature induced 
by the trap is the appearance of gaps in the linearization 
spectrum, which leads to the dependence of the largest 
unstable eigenvalue on C, as shown in the top panel of 
Fig. [51 Interestingly, while the envelope of this depen- 
dence traces a curve similar to that in Fig. ^ the gaps 
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FIG. 5: (Color online) The most unstable eigenvalue (top) for 
the ID lattice Skyrmion, and an example of the field configu- 
ration and its stability for C — 0.15 (bottom left and right) in 
the presence of the confining parabolic potential with Q — 0.1. 



lead to restabilization of the discrete Skyrmion in cer- 
tain intervals. For instance, the solution shown is stable 
for 0.113 < C < 0.115, where it is unstable without the 
trap. Hence, the parabolic potential offers (which is also 
true in the 2D and 3D cases) an additional mechanism 
of stabilization of the discrete Skyrmions. 

Conclusions. We have constructed discrete structures 



that emulate Skyrmions on nonlinear dynamical lattices. 
We have used the "hedgehog" ansatz for 2D (baby-) 
Skyrmions in continuum field models as a guide towards 
constructing the lattice solutions and extending them to 
both 3D and ID settings. Generally, the lattice Skyrmion 
is built herein as a vortex soliton in a complex field cou- 
pled to a dark radial soliton in a coupled real field. We 
have predicted stability limits of the lattice Skyrmions 
and identified the principal mechanism of their destabi- 
lization. The analytical prediction was verified by numer- 
ical computations, that corroborate the stability of the 
Skyrmions in the predicted region. We have also demon- 
strated that the evolution of unstable discrete Skyrmions 
leads to onset of lattice turbulence. A possibility of fur- 
ther stabilization of the Skyrmions by means of an exter- 
nal confining potential was highlighted too. 



It would be interesting to examine analogs of the struc- 
tures presented herein with higher values of the topolog- 
ical charge, as well as more complex patterns in 3D. Cre- 
ation of the lattice Skyrmions in binary BECs mixtures 
trapped in a deep optical lattice by means of available 
techniques appears to be within experimental reach. 
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